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Relaxation in time elapsed neuron network 
models in the weak connectivity regime 
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Abstract 

In order to describe the firing activity of a homogenous assembly 
of neurons, we consider time elapsed models, which give mathematical 
descriptions of the probability density of neurons structured by the 
distribution of times elapsed since the last discharge. Under general 
assumption on the firing rate and the delay distribution, we prove the 
uniqueness of the steady state and its nonlinear exponential stability in 
the weak connectivity regime. In other words, total asynchronous firing 
of neurons appears asymptotically in large time. The result generalizes 
some similar results obtained in [HKITl in the case without delay. Our 
approach uses the spectral analysis theory for semigroups in Banach 
spaces developed recently by the first author and collaborators. 
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1 Introduction 

In nervous systems, neuronal circuits carry out tasks of information transmis¬ 
sion and processing. Many neurons generate trains of stereotyped electrical 
pulses in response to incoming stimulations. Following each discharge, the 
neuron undergoes a period of refractoriness during which it is less responsive 
to inputs, before recovering its excitability Id- The main carrier of infor¬ 
mation is the discharge times or some statistics of the discharge times. In 
this work, we consider a simple neuronal model which neglects the mecha¬ 
nisms underlying spike generation and focusses on describing the neuronal 
dynamics in terms of discharge times. More precisely, we consider a model 
which has been introduced and studied in [3, |T6l [T7] and which describes the 
post-discharge recovery of neuronal membranes through an instantaneous fir¬ 
ing rate that depends on the time elapsed since the last discharge and the 
inputs by neurons. We refer to these papers for biologic motivation and dis¬ 
cussions. We also refer to n El 1221 ED] where these models (or similar ones) 
are obtained as a mean field limit of finite number of neuron network models. 

The neuronal network is described here by the density number of neurons 
/ = f{t,x) > 0 which at time t > 0 are in the state x > 0. The state of a 
neuron is a local time (or internal clock) which corresponds to the elapsed 
time since the last discharge. The dynamic of the neuron network is given 
by the following nonlinear time elapsed (or of age structured type) evolution 
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equation 


dtf = -d^f - a{x,em{t))f =: Cem(t)f, 
f{t,0)=p{t), f{0,x) = fo{x). 


(i.ia) 

(i.ib) 


Here a{x, e p) > 0 represents the firing rate of a neuron in the state x for 
a network activity /r > 0 and a network connectivity parameter e > 0. 
The function p{t) represents the total density of neurons which undergo a 
discharge at time t and is defined through 

pit) ■='P[f{t);m{t)], ( 1 . 2 ) 


where 


'P[9,p]='Pe[9,p]-= a{x,ep)g{x)dx. 


:i-3) 


The function m{t) represents the network activity at time t > 0 resulting 
from earlier discharges and is dehned by 


m{t) := / p{t - y)b{dy), 


:i-4) 


where the delay distribution 6 is a probability measure which takes into 
account the persistence of the electric activity in the network resulting from 
discharges (synaptic integration). In the sequel, we will consider the two 
following situations : 

• The case without delay, when b = and then m{t) = pit). 

• The case with delay, when 6 is a smooth function. 

We observe that in both cases, the solution / of the time elapsed equation 
fnyTj)-(iT7gj) satishes 


— / f{t,x)dx = f(t,d) 


a{x, emit)) fit, x)dx = 0. 


As a consequence, the total density number of neurons (also called mass in 
the sequel) is conserved and we can normalize that mass to be 1. In other 
words, we may always assume 


{fit, ■)) = ifo) = 1, Vt > 0, {g) := / g{x)dx. 
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A (normalized) steady state for the time elapsed evolution system of 
equations is a couple (F^, M^) of a density number of neurons 

Fe = F, {x) > 0 and a network activity > 0 such that 

0 = -d^F^ - a{x, e M^)F^ = C^m^Fs, (i-Sa) 

F,(0) = (F,) = 1. (1.5b) 

It is worth emphasizing that for a steady state the associated network activity 
and discharge activity are two equal constants because of the normalization 
of the delay distribution, i.e. {b) = 1. 

In equations f[T:Tj) -( iT7gj ) and ( [n^j ), the connectivity parameter e: > 0 
corresponds to the strength of the influence of the neuronal network activity 
on each neuron through the functions m{t) and f{t,x) respectively. In the 
limit case £ = 0, equation ffTTT| -( iTT^ is linear which means that each neuron 
evolves accordingly to its own dynamic. In the other hand, when e > 0, 
equation ffT7T| -( iTT^ is nonlinear and the dynamic of any given neuron is 
affected by the state (or the past states in the case of the model with delay) 
of all the other neurons through the global activity of the neuronal network. 
Finally, the weak connectivity regime, about which we are mainly concerned 
in the present paper, corresponds to a range of connectivity parameter e G 
(0,£o], ^0 > 0 small enough, such that the nonlinearity of equations IfTT]) - 
([1^ and (ji^ is not too strong. 

Our main purpose in this paper is to prove that solutions to the time 
elapsed evolution equation ffTTj) - ([r:^ converge to a stationary state under a 
weak connectivity assumption. Before stating that result, let us present the 
precise mathematical assumptions we will need on the bring rate a and on 
the delay distribution b. 

We make the physically reasonable assumption 

dxtt >0, a' = d^a > 0, (1.6) 

0 < Oo := lim a(a;, 0) < lim a{x,fi) =: ai < 00, (1.7) 

X,{1^00 

as well as the smoothness assumption 

a e W 2 ’'^(M^). ( 1 . 8 ) 

In the delay case, we assume that b{dy) = b{y) dy satisfies the exponential 
bound and smoothness condition 

poo 

35 > 0 , / {b{y)+ \b'{y)\)dy < 00 . ( 1 . 9 ) 

Jo 
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We begin by stating our main result about the stationary problem ( [TTgj ). 

Theorem 1.1. Assume fli.6l) -( iTT7j )- fli.8l) . For any e: > 0, there exists at least 
one solution {F^, M^) G x M_|_ to the stationary problem such 

that 

an . 

0 < F,(x) < Ce-2^, \Fl{x)\ < C e--^^, Vx > 0, (i.io) 

for a constant C G (0, cxo). Moreover, there exists Eq > 0, small enough, such 
that the above solution is unique for any e G [0,£o)- 

For a given initial datum 0 < /o G L^(R+), we say that a function / is a 
weak (positive and mass conserving) solution to (fT7n) -( iTT^ if 

0 < / G C'([0,oo);Li(M)), (/(f)) = (/o), Vf > 0, 

and / satisfies (fTTTl) in the distributional sense P'([0, cxd) x [0, oo)) for some 
functions m,p E ^([O, oo)) which fulfilled the constrains (fiT^) and ( [TT^ . 

Under the above assumptions, existence and uniqueness of weak solutions 
have been established in |2^ Theorem 1.1], The main concern of the present 
work is the following long-time asymptotic result on the solutions. 

Theorem 1.2. We assume that the firing rate a satisfies dUnD, (iiTl) and 
(li.Sh . ITe also assume that the delay distribution b satisfies b = do or ■ 
There exists Eq > 0, small enough, and there exist some constants a < 0, 
U > 1 and T] > 0 such that for any connectivity parameter e G (0,eo) 
any initial datum 0 < fo E with mass 1 and such that ||/o — < rj/E, 

the unique solution f to the evolution equation IfTTlj) -( fTT^ satisfies 

||/(f,.)-F,|Ui <C'e“^ Vf >0. 

In other words, in that weak connectivity regime, we prove that the total 
asynchronous firing of neurons appears exponentially fast in the large time 
asymptotic. Theorem 11.21 extends to firing rates a satisfying fli.6|) - fli.8p some 
similar exponential stability results obtained in [T6l ITT] in the case without 
delay and for a step function firing rate a given by 

a(x,/i) = 1„><^(^), cr,a"^ G 1U^’°°(R+), a'< 0. (i.ii) 

It is worth mentioning that the above firing rate does not fall in the class of 
rates considered in the present paper because condition (li.8p is not met. On 
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the other hand, we are able to tackle the case without and with delay in the 
same time, what it was not the case in In the delay case, stability 

results were established in [16], but not exponential stability. 

Our proof follows a strategy of “perturbation of semigroup” initiated in 
[ 12 ] for studying long time convergence to the equilibrium for the homoge¬ 
neous inelastic Boltzmann equation and used recently in [13] for a neuron 
network equation based on a brownian (hypoelliptic) perturbation of the 
well-known FitzHugh-Nagumo dynamic. More precisely, we introduce the lin¬ 
earized equation for the variation functions {g, n, q) = (/, m, p) — {F^, M^, M^) 
around a stationary state {F^, M^, M^), which writes 


dtg = -d^g - a{x,eMs)g - n(t) e{d^a){x,eMe)Fe, ( 1 . 12 a) 

^(f, 0 )=g(f), ^( 0 ,a;) = ^o(a;), ( 1 . 12 b) 


with 



a(a:, £ Mi;)g dx + n(t) e 


{df,a){x,eMs)F^dx 


(1-13) 


and 


n{t) := / q{t-y)b{dy). 


( 1 . 14 ) 


^'O 

We associate to that linear evolution equation a generator (which acts on 
an appropriate space to be specified in the two cases without and with delay) 
and its semigroup Sa^. It turns out that we may split the operator as 


+ ^61 


for some a-hypodissipative operator a < 0 , and some bounded and 
power regular operator Ae as defined in [211 [HI [HI [TUI lH]- particular, 
adapted versions of the Spectral Mapping Theorem in [HI [TOl [TT] and the 
Weyl’s Theorem in [211 El [HI [TUI E] imply that the semigroup S'a^ has a fi¬ 
nite dimensional dominant part. Moreover, in the limit case when e: = 0, the 
term n(t) disappears from equation In. 121 ) and the resulting semigroup S'a^ 
becomes positive. That allows us to use the Krein-Rutman Theorem estab¬ 
lished in [HI [ini E] in order to get that the stationary state (Fq, Mq, Mq) is 
unique and exponentially stable. Using next a perturbative argument devel¬ 
oped in [121 [231 El E], we get that the unique stationary state (F^, M^, M^) 
is also exponentially stable in the weak connectivity regime. We conclude the 
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proof of Theorem II. 2l bv a somewhat classical nonlinear exponential stability 
argument. 

The same strategy applies to the case without delay and with delay. In 
both cases, the boundary condition in the age structure equation is treated 
as a source term and, in the delay case, the delay equation ( [TT^ is replaced 
by a simple age equation on an auxiliary function, so that the resulting 
linearized equation writes as an autonomous system of two PDEs and falls in 
the classical framework of linear evolution equations generating a semigroup. 

Our approach is thus quite different from the usual way to deal with delay 
equations, as introduced by I. Fredholm [6] and V. Voltera [25], which consists 
in using the specific framework of “fading memory space", which goes back 
at least to Coleman & Mizel [I], or the theory of “abstract algebraic-delay 
differential systems" developed by O. Diekmann and co-authors |3]. 

Our approach is also different from the previous works usiiniiiH] where 
the asymptotic stability analysis were performed by taking advantage of the 
step function structure fii.iii) of the firing rate. That one makes possible 
to explicitly exhibit a suitable norm (related to the Wi Monge-Kantorovich- 
Wasserstein optimal transport distance) such that some related linear age 
structure operator is dissipative. The present method is based on a more 
abstract approach but in the other hand it is somewhat more flexible because 
if does note require to explicitly exhibit a norm for which the underlying 
linear(ized) operators are dissipative. In particular, we hope that our strategy 
can be adapted to the large connectivity regime, to the step function firing 
rate fii.iii) as well as to models including fragmentation term to describe 
neuronal networks with adaptation and fatigue, and thus generalize to the 
case with a delay term all the stability results established in [TGl IT71 118] in 
the case without delay. 

Let us end the introduction by describing the plan of the paper. In 
Section 121 we introduce the strategy, we prove the stationary state result and 
we establish Theorem ll.2l in the case without delay. In Section[3l we establish 
Theorem 11.21 in the case with delay. As we mentioned above, the strategy of 
proof for the case without and with delay is rather the same. For pedagogical 
reason, we start presenting the method on the simplest “without delay case " 
in Section [2l where we prove the stationary problem result Theorem 11.11 as 
well as Theorem ll.2l in that case. Next, in Section[3l we only explain how the 
proof must be modified in order to treat the more complicated “with delay 
case " and thus establish Theorem 11.21 in all generality. 
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2 Case without delay 

The present section is devoted to the proof of our main result Theorem 11.21 
in the without delay case. 

2.1 The stationary problem 

We first deal with the stationary problem and we prove the existence of 
steady state as well as its uniqueness in the small connectivity regime. 

Proof of Theorem M . 1{ Step 1. We prove the existence of a solution. We set 

nx 

A{x,m) := / a{y,m)dy, Vx, m>0. 

Jo 

For any m > 0, we can solve the equation ( ii. 5 aj ), by writing 

( 2 . 1 ) 

where > 0 is chosen in order that satisfies the mass normalized 
condition, namely 

poo 

T-i = / 

Jo 

In order to conclude the existence of a solution, we just have to find a real 
number m = such that m = F^^rn{d) = Tm- Equivalently, we need to find 
Me > 0 such that 

$(£,MJ = 1 , ( 2 . 2 ) 

where 

poo 

<F(£,m) = rnT~^ :=m 

Jo 

From the assumption (Eli) of a, there exists Xq G [ 0 , 00 ) such that 
a(x, p) > Yy ^ Xq, p > 0, and therefore 

^(x — xo)+ < A(x, p) < aix, Vx > 0, V/x > 0. ( 2 . 3 ) 
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Thanks to the Lebesgue dominated convergence theorem, we dednce that 
<h(£, •) is a continnons fnnction. Becanse $(^,0) = 0, <h(e:, oo) = oo and 
thanks to the intermediate valne theorem, we conclnde to the existence of at 
least one real nnmber G (0, oo) snch that holds. Estimates fii.ioi) 
immediately follow from the identity feTT]) and the estimate (j3^. 

Step 2. We prove the uniqueness of the solution in the weak connectivity 
regime. Obvionsly, there exists a nniqne Mq := (f^ G (0, oo) 

snch that <h(0,Mo) = 1. Moreover, we compnte 

(9 /■°° dA 

— —<h(e,m) = / — me- — (x,em))dx, 

UTTl Jq UTTl 

which is continnons as a fnnction of the two variables becanse of fli.81) . We 
then easily obtain that <h G Since moreover 

8 /-oo 

—^{e,m)\e=o = j > 0 , 

the implicit fnnction theorem implies that there exists eo > 0 , small enongh, 
snch that the eqnation fe:^) has a nniqne solution for any e G [0,eo)- Q 

Remark 2.1. In the above proof, we do not need fli.Sp but only the weaker 
smoothness assumption that A and d^A are continuous. 

2.2 Linearized equation and structure of the spectrum 

To go one step further, we introduce the linearized equation around the 
stationary solution On the variation {g,n), the linearized equation 

writes 


dtg + 9x5' + asg + a'^Feulf) = 0, 

poo 

g{t, 0) = n{t) = / {a^g + a'^F^n{t)) dx, 5 (( 0 , x) = go{x), 

Jo 

with Oe := a{x,eM^), a'^ := £ {d^a){x,eM^). Since there exists Eq > 0, small 
enough, such that 


Ve G (0,eo) 


K : = 


a^Fgdo: < 1, 
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we may define 


r 

■M.e[g\ := (1 - / 

^0 


'OO 


a^g dx, 


(2.4) 


and the linearized equation is then equivalent to 


dtg + dxg + a^g + Me[g{t ,.)] = 0 , 

= Me[g{t,.)], g{Q,x) = go{x). 


(2.5) 

(2.6) 


To the above linear evolution equation, one can classically associate a 
semigroup Si^it) acting on X := with generator 


Leg ■= -dxg - a^g - Me[g] 


and domain 


D{L,) := {g G g{Q) = MM)- 


Here we also use another approach by considering the boundary term as 
a source term, and then rewriting the equation as 


dtg = Aeg := -dxg - aeg - a'^FeMe[g] + 5 x=oMe[g], (2.7) 

acting on the space of bounded Radon measures 

A” := = {g e (C'o(M))'; supp^f C M+}, 

endowed with the weak * topology cr(M^, Co). Here and below, Co(/) denotes 


the space of continuous functions on the closed interval J = M or J = M+ 
which goes to 0 at inhnity. We then denote by SA^{t) the semigroup on A 
generated by A^. 

It is worth emphasizing that for any gQ G X, the function g{t) = Sl^ {t)go £ 
C([ 0 , cxo);X) is a weak solution to equations ( p75j )- fl2.6p . and more precisely 
clearly satishes 



for any (p G C^([0, oo); Co(M_|_)) fl Cc([0, cxo); Cq(M+)). Here, we have dehned 
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for any G C'q(M+), the space of functions which goes to 0 at infinity 
as well as their first derivative. As a consequence, the semigroup S'a^ being 
defined by duality from the semigroup we have S'aJv = Sl^. 

For a generator L, we denote by S(L) its spectrum and hy Sl the as¬ 
sociated semigroup. We refer to the classical textbooks p ca m for an 
introduction to the spectral analysis of operators and the semigroup theory. 
Our next result deals with the structure of the spectral set ^(Ag) of A^ and 
the splitting structure of the associated semigroup S\^. 

Theorem 2.2. Assume fli. 6 p -( iT? 7 j )- (li. 8 p and define a* := —ao/2 < 0. The 
operator A^ is the generator of a weakly * continuous semigroup Sa^ acting 
on X endowed with the weak * topology cr(M^,Oo). Moreover, there exists a 
finite rank projector nAj,a* which commutes with Sa^, an integer j > 0 and 
some complex numbers 

G := {zeC, ^ez> a*}, 

such that on Ei := YIas^wA^ the restricted operator satisfies 

S(A,|sJnA„* = (2.8) 

(with the convention S(A£|eJ O A^* = 0 when j = Oj and for any a > a* 
there exists a constant Ca such that the remainder semigroup satisfies 

||^A,(/ - nA„a‘)IUw < Vt > 0. ( 2 . 9 ) 

The proof of the result is a direct consequence of the fact that the operator 
Ae splits as K,, = Ae + where and are defined on X by 

Aeg := yeMe[g], 7 ^ := (5o - a'F^, ( 2 . 10 ) 

Beg :=-d,,g - Oeg, ( 2 . 11 ) 

for which we can adapt Weyl’s Theorem of IMlElElCnilll] and the Spectral 
Mapping Theorem of [HCniE]- The picture may seem not that simple, 
because Ag does not generate a strongly continuous semigroup on X and 
then apparently does not fit the framework developed in [TTl [TO]. However, 
we may probably circumvent that issue in the following ways: 

- We may observe that the strong continuity is little used in [TTl [10] and 
then the results stated therein extend to weakly * continuous semigroups. In 
other words, on the finite dimensional eigenspace associated to the principal 
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part of the spectrum fl 2 . 8 l) continuity and weak * continuity are equivalent, 
while in the remainder part ( p:gj ) we just use a decay bound which does not 
require the strong continuity property, see [T^ Chapter 1 ] as well as [261 [TT] 
where such a weak * continuous framework is also discussed. 

- We may apply the theory developed in naiin] to the adjoint operator 
A* and the associated semigroup acting on Co(M) and then deduce the result 
by duality. 

- We may probably reduce the space X by “sun duality", see jH Chap¬ 
ter II.2. 6 ], apply the theory developed in [HI [10] to the resulting strongly 
continuous semigroup and then conclude by a density argument. 

We rather follow another strategy, which has the advantage to be more 
self-contained and pedagogical, by adapting to our context the proofs from 
da [ID]. It is worth emphasizing that we introduce the two spaces X (Z X 
because the splitting has some clearer meaning in the space 

X while the Banach structure of X is used when we establish the adequate 
version of Weyl’s theorem. 

We recall the definition of hypodissipativity introduced in [ 8 ]. We say 
that the closed operator L on a Banach space E with dense domain D{L) is 
a-hypodissipative if there exists an equivalent norm ||| ■ ||| on E such that 

WfeD{L), A|||.|||(/) 3fJe(^,(L-a)/)<0, 

where, for any f G E, the associated dual set F|||.|||(/) C E' is defined by 

A|||.|||(/) := {V. e E'; (^,/) = 111/111^=111^1111.}. 

We also recall that the generator L of a semigroup of bounded operators is 
a-hypodissipative if, and only if, there exists a constant M > 1 such that 
the associated semigroup Sl on E satisfies the growth estimate 

||5'L(t)||^(^;) < Me“/ Vf>0, 

where 3S{E) denotes the space of linear and bounded operators on E. We 
will sometime abuse by saying that Sl is a-hypodissipative when it satisfies 
the above growth estimate. We refer to da El El n] for details. 

We start with the properties of the two auxiliary operators. 

Lemma 2.3. Assume that a satisfies conditions fli.6p - fli.8p . For any e > 0, 
the operators and satisfy : 
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(i) A, G 

(ii) Sbs is a*-hypodissipative in both X and X. 

(Hi) The family of operators Sb^ * A^Sb^ satisfies 

\\{SB,*AeSBmyix,Y) < Cae^\ Va > a*, 

for some some constant Ca G (0, oo) and with Y := nLj(M+). 

Here i?y(M+) stands for the space of functions with bounded variation and 
LJ(M_i_) stands for the weighted Lebesgue space associated to the weight func¬ 
tion XI—)- (x). 

Proof. In order to shorten notation we skip the e dependency, we write a(x) = 
a{x,eM^), A{x) = A{x,eM^), S'a = S'a^, S_a = Sb = Sb^ and so on. 

Step 1. Proof of (i). We have A G df) from the fact that 

G M) because ||a||wi.°° < cxo as assumed in fli.Sp . 

Step 2. Proof of (ii). We use the same notation B and Sb for the generator 
and its associated semigroup dehned in the spaces X and X. We write Sb 
with the explicit formula 

SB{t)g{x) = - t)l^_t>o =: S{t, x). ( 2 . 12 ) 

From the inequality 0 ( 2 :, 0) > (3ao/4)(l — lo<z<xi) for some Xi G (0, 00 ) 
coming from ( [TTTj ), we deduce A{x — t) — A{x) < 3aoXi/4 — (3ao/4)f for any 
X > f > 0, next 

^A{x-t)-A{x) < Vx > f > 0, ( 2 . 13 ) 

3ar\x-\ 

with C := e~^ > 0 , /3 := —ao/4 < 0, and hnally 

\\SBit)g\\x<Ce^^^\\g\\x, Vf > 0, ^geX. 

We conclude by observing that 3/3 < a* and that the same estimate holds 
with X replaced by X thanks to a (weakly *) density argument. 

Step 3. Proof of (Hi). For g G X and with the notation of Step (ii), we have 

ASB(t)g = '7eN(t), 
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where 7 ^ is defined in fi 2 .ioj) and 


POO 

N{t) := Me[S{t, •)] = (1 - k)~^ / dx. 

Jo 

It is worth noticing that N G C6(M+), because {x,t) 1—)■ is 

a bounded and continuous function on the set T := {{x,t)] x > t > 0 }. 
Moreover, as in Step (ii), we have 

POO 

\N{t)\<Cai e^^^\g{x-t)\l^-t>odx <Caie^^^\\g\\x, (2.14) 

Jo 

for any f > 0 . We deduce 

{SB*ASB){t)g{x) = [ {SB{s)%){x)N{t - s)ds 

Jo 

^A(.-s)-Ai.) _ s)l,_,^,ds (2.15) 

= {Ue* Nt){x), (2.16) 

with := ''fse^ and the classical notation W(s) = N{t — s). Starting from 
identity ( 12.151 ) denoting e-i3{x) := e~J^, we have 

|(S'B*MS'B)(f)5((a;)|e_/}(a;) < f C {'jeix - s)\\N{t - s)\ ds 

Jo 

<C'e“‘*||(| 7 s|e_/,)*(|iV|e_^)||c, 

< Ce“**|| 7 £e_^||A' \\Ne-fs\\c,- 

Using the definition fi2.ioj) and the estimates f|i.ioj) and ( 12.141 ), deduce 
||( 5 b < C"e“** ll^ll;^, Vf> 0 , (2.17) 

for some constant C G (0, 00). 

Next, differentiating the functions in both sides of identity fl2.i6p . we get 

dx [(S'e * ASb){ t)g] (x) = -a(x)e~^‘^^^ (z/^ * W)(x) - (z/^ * N')(x) 

I^eix - t)N(0)lx-t>o + i'e{x)N{t), 



14 














with N^{s) = N'{t — s). In order to estimate the second term, we compute 

roo 

N\t) = (1 - dx 

Jo 

poo 

= [a'{x) - -t)l^_t>odx. 

Jo 

Using the inequality ( i 2 .i 3 j ), we deduce 

poo 

|^^(^)| < C' {|a'(x)| + a(a:)^} — t)l3;_i>o dxdt 

< e^*'\\g\U- 

As a consequence, using again ( | 2 .i 3 j ), we get 

pt poo 

NOWx ^ J J s)\\N'{t-s)\dxds 

Treating in a similar way the two other terms involved in the expression of 
dx[{St 3 * ASB)it)g], we finally obtain 

m{SB*ASB)mh<e'^^^\\9h. 

Together with that concludes the proof of (hi). □ 

Proof of Theorem \2A The main idea is to apply or adapt the versions of 
Weyl’s Theorem m Theorem 3.1], (lOl Theorem E.3.1] and of the Spectral 
Mapping Theorem m Theorem 2.1], [TOl Theorem E.2.1] (see also the vari¬ 
ant results in m)- We start collecting the three key properties satisfied by 
the operators As and Bs involved in the splitting of 

We denote by the abstract Sobolev space associated to Ls for C ^ 
see e.g. [H Section II. 5], so that in particular C IU‘’’^(M+), with equality 
when C < 0. We recall that for the generator L of a semigroup Sl, we define 
the resolvent set p{L) by 

p{L) := { 2 ; G C; L — z is a bijection} = C\S(L), 

as well as the resolvent (operator) Rl{z) := (L — z)~^ for any z G p{L). We 
finally define the half plane := {z G C; ^ez > a} for any a G M. 
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(Al) For any a > a* and £ G N, there exists a positive constant Ca/ such 
that the following growth estimate holds 

11% * < Ca,,e^\ yt > 0. (2.18) 

It is obvious that (Al) is true for i = 0 from Lemma [2.31 — and for 
^ = 1 from (\2.iY^ - We then deduce (Al) for any i >2 and a > a* by 
writing % * (Ae%)^*^^ = [% * (Ae%)] * (Ae%)(*^“^) and by using 
that {Be — a*) is hypodissipative in X and Ae G ^(A). 

(A2) There holds As G e^(A_i,A^) for any C, G (—1,0) and the family of 
operators (%Ae)*'*^^(f) satishes the estimate 

||(%Ae)^*^^(f)||^(x_i,x) < \ Vf > 0, (2.19) 

for a positive constant The hrst claim is a consequence of the 

continuous embedding X C and of Lemma [2.31 — (z). The second 
estimate is obtained by putting together the properties {%) and {ii) in 
Lemma 12.31 

(A3) The family of operators {R! 3 ^AeY{z) satishes the compactness estimate 
\\{RBe{^)^ef‘\\.'S{X,Y) — ^ 2 ,ai V 2 ; G A^, \/ a > a*, 

for some positive constant Observing that i—)■ —{Rs^{z)AeY is 

nothing but the Laplace transform of the function t i—)■ (S'g^Ae)^*^^(f), 
that clearly holds true thanks to Lemma [2.31 

We now briehy explain how the proof goes on, referring to [Til Uni Cl] for 
more details and developments. In order to shorten notation we skip again 
the £ dependency. 

Step 1. On the one hand, because of ( |2.l9j ) and the fact that Rs is nothing 
but the opposite of the Laplace transform of Sjg, we have 

\\{Rb{z)AY\\x_i^x < c, WzeAa. ( 2 . 20 ) 

Recalling that the negative abstract Sobolev norm is dehned by 

||5'||x_i := \\RB{a)g\\x, 
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and using the dissipativity property of B, we immediately have 

Vz e Aa, \/g E X, \\RQ{z)g\\x_i < C ||5'||x_i- 

Moreover, using the resolvent identity 

Rb{o)Rb{z) = z ^ {RB{a)BRi3{z) — RBia)), 

= z~^ (RBiz) - RBia) + aRB{z)RB{a)), 


we also have 

||-Rb(^)5'IU_i < -p-i {||-Rb(^) 5'IU + ||-Rn(a)5'l|v + |a| ||-RB(^)-RB(a)5'l|x} 

< mllsllx. 

for any g E X and E A^. Finally, using both estimates and an interpolation 
argument, we deduce 

VzeA„, \\Rb{z)\\x^^x.i < VC e (- 1 , 0 ). ( 2 . 21 ) 

We conclude with 

^ZE Aa, \\{RB{,z)Af\\x^x < ( 2 - 22 ) 

by writing 

{RB{z)Af = {RB{z)AfRB{z)A 

and by using estimates (i) in Lenima [273] with C = —1/4, 112 . 21 ) with C = —1/4 
and 112 . 20 ) . 

Step 2. Here, we follow mi Proof of Theorem 2.1] and |24) by taking 
advantage of the analysis of degenerate-meromorphic functions performed in 
|21) . Iterating the factorization identity = Rb ~ Rb^Ra, obtain 

{I-V)Ra=U in ^{X) ( 2 . 23 ) 

with 

U = R]b — RbARb ~\~ {RbA)^Rb, V •= — {Rb-^)^- 
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Because of fi 2 . 22 j) . fb.iSp and of the identity 

poo 

RBiARsYiz) = (- 1 )^+^ / e-^^Ss * dt e ^(X) 

Jo 

for any 2 ; G A^, we see that all the terms involved in (gT^jj) act in ^(X), so 
that we also have 

{I-V)Rl=U in SS{X). ( 2 . 24 ) 

Because of fl 2 . 22 i) . we see that / — V is inversible on A(j\i?(0, M) for M > 0 
large enough and ( 12 . 24 ! ) tells us that Rl = (/ — is uniformly bounded 

in Aa\B{0, M), in particular S(L) C A^ U 5(0, M). 

On the other hand, because of (A3), there holds 

||V(2;)||x^y < C, 

with compact embedding Y G X. As a consequence of the fact that V is 
the opposite of the Laplace transform of the nice function R+ —)■ A, t 1 —)■ 
(SBAG>(t), we also deduce that V is holomorphic in ^{X) on A^. We 
define := I — V{z). Because d' is an holomorphic function on Aa, 'L(^) 
is invertible for (some) z G Aa\B{0, M) and RV{z) C Y for any 2 ; G A^, 
[ 21 I Corollary II] asserts that the function ^ \['(z)“^ is a degenerate- 

meromorphic operators valued function, that is d' is holomorphic on Aa\D, 
D C Aa is discrete, any point ^ G 5 is an isolate pole and the coefficients of 
the principal part in the Laurent series are finite rank operators. Next, from 
[ 21 ] and the identity 

Rl = 

we deduce that Rl is also degenerate-meromorphic on A^, in particular 
S(L) n Aa is a finite set of discrete eigenvalues, which means isolated eigen¬ 
values associated to an algebraically finite eigenspace. That is nothing but 
fb.Sp . We define If := / — IlA.a- 
Step 3. We claim that for any n > 2 , there holds 

SLm = { 5 B(t) + ... + (^BA)*("-'U 5 B(t)}n (2.25) 

+(- 1 )"+'^ / (RBAnz) RL{z)Ue^^ dz 

in B{X) where tc:= {c + iy, y G M} is the complex line of abscissa c G M. In 
order to prove that identity, we iterate the Duhamel formula 

5A(t) = 5B(f) + 5BA*5A(f), 
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in and for any / G X and n>2. We then obtain 


SLmf = {Ss{t) + ... + * Ss{t)}Uf + gn{t), 


with 

gnit) := * SLmf- 

Because t i—)■ g'^{t)e~^^ = Bgn{t) e~’^^ + Agn-i{t) G /.^(R+jX) for 6 > 0 
large enough from (A1)-(A2) and the fact that 1 1-)- S'L(t)n/e“^* G L°°{X), 
the following inverse Laplace transform formula holds 

= [ (RBAnz) RL{z)uf dz, vf>o. 

2vr 

From the preceding steps, the function under the integral sign is holomorphic 
on Aq* and we may move the line of integration from to ta- Observing 
that every term in the resulting identity belongs to X, we conclude to (\2.2^^ . 
More precisely and importantly, taking n = 6, using (i2.22j) and the fact that 
\\Rl{z)I\\\^(^x) is uniformly bounded on from Step 2, we have 

at r 

\\ 9 n{t)\\x < — / \\iRBmi^)\\.^ix)\\RLiz)'A\\^(x)dz\\f\\x < e°-^\\f\\x, 

ta 

for any f > 0 . We conclude to using that the other terms involved in 
(\2.2^^ are similarly bounded thanks to (Al) and using a weakly * density 
argument. □ 

2.3 The vanishing connectivity regime 

When the network connectivity parameter vanishes, e = 0 , the linearized 
time elapsed operator simplihes into 

Aog = -d^g - a(x, 0 )g + d^=QMo[gf (2.26) 

where M.Q[g] = a(x, 0)5f(x)da:. The associated semigroup is then positive 
and the following version of the Krein-Rutman theorem holds. 

Theorem 2.4. There exist some constants Oq < 0 and C > 0 such that 
S(Ao) n Aqq = { 0 } and for any g^ G X, (go) = 0 , there holds 

\\SAo{t)go\\x < Ce'^^^WgoWx, Vt > 0 . (2.27) 
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We denote 


A’+ := {g e M\R+)-, g > 0}, 

the space of bounded and nonnegative Radon measures. 

We start with two elementary auxiliary results. 

Lemma 2.5. S'aq is positive: SAo(t)g ^ for any g G and any t > 0. 
Proof. We introduce a dual problem of fl 2 . 26 l) defined on the space C'o(M) by 


dt(p = A(p = Bip + Aif, <p(0, •) = (po, 
where the operators A and B are defined by 

Bip = dx(p — a{x,0)ip, yl,<p = a(a;, 0)<p(0). 
A solution ip to fl 2 . 28 p then satisfies 

ip{t) = S^{t)ipo + {S^ * Aip){t). 


( 2 . 28 ) 


Let us fix cpQ G Co(M) such that (po < 0 and let us prove that (p{t) < 0 
for any f > 0. We obviously have that is a positive operator and it is a 
contraction in Co(M). Taking the positive part in fl 2 . 28 p . we get 

P+it) < S^{t)(po+ + (Sg * A(p+){t) 


so that 


<ai S^{t - s)p:+{0)ds, 


\\p:+{t)\\L^<C ||v 7 +(s)||l.x= ds. 


From Gronwall’s lemma, we deduce that (p+(t) = 0 for any t > 0 and then 
(p < 0. We conclude by observing that Saq is the dual of Fy. □ 

Lemma 2.6. —Aq satisfies the following version of the strong maximum 
principle: for any given g G A’_|_ and /i G M, there holds 

g G L*(Ao) \ {0} and (—Aq + p)g > 0 imply g > d. 
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Proof. Suppose that there holds (—Aq + ^)g > 0 for g satisfying the above 
conditions. It is only necessary to prove that g does not vanish in M_|_. Since 
g ^ 0, there exists x* G M+ such that g{x*) > 0. Rewriting the assumption 
as 

poo 

dxg + {a{x, 0) + g)g > 5 x=q / a{x, Q)gdx, 

Jo 

we observe that 

^ g^(a., 0 )+Mx 0 ) + /i)^) > 0 . ( 2 . 29 ) 

(i) For X G {x*, 00 ), we have 

^A{x,o)+^,Xg > > 0 . 

(ii) For X G (0,a:*), by integrating the same equation on (0,a:), we obtain 


^A{x,0)+t.xg > / 


a{y,0)g{y)dydu + g{0) 


> 


a{y,d)9{y)dy. 


From the positivity assumption (fTTyj) on a and step (i), we have 


a{y,d)g{y)dy > y 


' max{a:o3^*} 


g{y)dy > 0. 


Therefore, g does not vanish on (0, 00 ). □ 

Proof of Theorem 2-4■ First, we know from Theorem 11.11 that there exists 
at least one nonnegative and non-vanishing solution Fq to the eigenvalue 
problem AqFo = 0 and the associated dual eigenvector is -0 = 1- Next, we 
observe that, dehning the sign/ operator for / G T*(Aq) by 

[(sign/)*/>](a;) := [f{x)f’{x) + f{x)i>{x)], G C'o( 

we have, for any -0 G C'o(M)+, 

3?e((sign/)ATo[/] ,/>) = 3fJe(ATo[/], (sign/)» 

3?e/(0) 


= 3?e 


aof dx 


^( 0 ) 


< 


l/( 0 )| 

ao\f\dx^p{0) = (Ado[|/|] ,-0), 
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which is nothing but the complex Kato’s inequality 

V/gD(Ao), 3fJe(sign/) Ao/< Ao|/|. ( 2 . 30 ) 

We also observe that -D(Aq) C (^^(OjCX)), and, as a consequence, g G -D(Aq) 
and > 0 implies g > Q 01 g < ^. We then may use exactly the same 
argument as in [TH IT 0 [ Proof of Theorem 5.3] (see also |TT|): 

- Kato’s inequality ( 12 . 301 ) and the strong maximum principle imply that 
the eigenvalue A = 0 is simple and the associated eigenspace is Vect(Fo); 

- together with the fact that Saq is a positive semigroup, one deduces 
that A = 0 is the only eigenvalue with nonnegative real part. 

We conclude to the spectral gap estimate ( 12 . 27 ] ) some Oq G (a*,0) 
with the help of Theorem 12.21 □ 

2.4 Weak connectivity regime - exponential stability of 
the linearized equation 

We extend the exponential stability property which holds for a vanishing 
connectivity to the weak connectivity regime thanks to a perturbation argu¬ 
ment. 

Theorem 2.7. There exist some constants £0 > 0, a < 0 and C* > 0 such 
that for any e G [0, £ 0 ] there hold ^(Ag) n Aq, = {0} and 

||*5Ae(^)^olU < C'e"* ll^olU, Vf>0, ( 2 . 31 ) 

for any go G T, {go) = 0. 

The proof follows the stability theory for semigroups developed in Kato’s 
book |9] and revisited in [121 ISSl IHl ITS] . 

Proof of Theorem \2. 71 With the definitions (g^j), 112 . 101 ) and Ij^.iij) of 
Ae and Sg, we have 


{Be - Bo)g = (a(x, 0) - a(x, eMe))g 


and 


{Ae- Ao)g = {Me[g] - Mo[g]) 5o - e{df,a){x,eMe) Fe Me[g]. 


22 














Together with the smoothness assumption fli.Sj) . we deduce 

\\Be-Bo\\^{x) + \\As-Ao\\ss{x)<Ce, Ve:>0. ( 2 . 32 ) 

We then argue similarly as in the proof of |231 Theorem 2.15] (see also [HI [121 
[mils]) and therefore just sketch the proof. We now dehne 

K,{z) := {RsMAefRAoi^KAs - Ao) e 

and we take any a G (ao,0), recalling that a* < ao < 0. We deduce from 
( 12 . 32 ] ) and the estimates (i) and (ii) in Lemma 12.31 that for some rj < |q;|, 
£0 > 0, small enough, C > 0, and for any G Aa\B{0,r]), e G [ 0 , 60)5 we 
have 

\\Keiz)y(x) <Ce <Ceo<l, ( 2 . 33 ) 

and thus (1 — Ks{z))~^ is well dehned in ^{X). Using the elementary iden¬ 
tities 

RAs = RBs — RSeAeRBe + {RBeAef Rke ='■ + {RBsAsY Rke, ( 2 . 34 ) 

and 

Rke = Rao + -RAo(Ae — Ao)i?A^, 

we get 

(/ - =U, + {RbAs^Rao- 

Because both terms are well dehned in ^{X) from Lemma [2.31 and the link 
between resolvent and semigroup, we also have 

{I-K,)RL,=U, + {RB,AefRLo, 
from what we deduce 


Rl, = (J - K,)-\U, + {RbAsYRao)- 


Since the RHS expression is clearly uniformly bounded in ^^{X) on the com¬ 
plex set AQ,\i?( 0 ,? 7 ) for any e G [ 0 , 60)5 we have ^(Lg) fl C 

Using the dehnition of the eigenprojector Llo on the eigenspace associated 
to the spectral values of Ao lying in B{0, rj) by mean of Dunford integral (see 
[9l Section III. 6 .4] or [H [TT] ) and the analogous of ([ 2 . 34 ] ) in AS{X), we have 


Un = 


2ti 

i 


'\z\=V 


’\z\=r) 


Rlo (z) dz 

) 

{RbqAoYRlq dz, 
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by using that the contribution of holomorphic functions vanish. In a similar 
way, we have 



For g G X, we next compute 




and the estimates and As a consequence, we deduce 

lin^ — no||i^(js:) < 1- 

for any e: G (0,eo); up to take a smaller real number £o > 0. From the 
classical result m Section 1.4.6] (or more explicitly [231 Lemma 2.18]), we 
deduce that there exists G Aq, such that 


^(Ag) n Aq, = {^e}, is a simple eigenvalue. 


for any £ G [0, £o]- We conclude by observing that Ce = 0 
A*1 = 0 (which is nothing but the mass conservation). 


0 because 1 G X' and 
□ 


2.5 Weak connectivity regime - nonlinear exponential 
stability 

Now, we focus on the nonlinear exponential stability of the solution to the 
evolution equation ffTTTi) -( [TT^ in the case without delay. We start with an 
auxiliary result. We define the function 4) : L^(M_|_) x M —)• R by 
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We denote by Wi the optimal transportation Monge-Kantorovich-Wasserstein 
distance on the probability measures set P(M+) associated to the distance 
d(x, y) = \x — y\ f\ 1 , or equivalently dehned by 



Lemma 2.8. Assume fli.Sp . There exists eo > 0 for any e G (0,eo) 
there exists a function : P(M) —)• M which is Lipschitz continuous for the 
weak topology of probability measures and such that p = (pe[g] is the unique 
solution to the equation 


/r G M+, ^{g,p)=Q. 

Proof of Lemma \2.8l Step 1. Existence. For any g ^ P (M) we have 
^{g, 0) > 0, and for any q G P(M) and /i > 0, we have 


^{g,T) < l|a|U- - T, 


so that ^{g,p) < 0 for /i > ||a||2,oo. By the intermediate value theorem and 
the continuity property of <h, for any hxed g G P(M+) and £ > 0, there exists 
at least one solution p G (0, ||a||Loo] to the equation ^{g,p) = 0. 

Step 2. Uniqueness and Lipschitz continuity. Fix f,g G P(M+) and consider 
/r, z/ G M+ such that 


<F{f,fi) = <h{g,n) = 0. 


We have 



with 



and 



We then obtain 


l/i - z/| (1 - e\\d^a\\L^) < \\a{-,en)\\wh^ Wi{g,f), ( 2 . 35 ) 
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and we may fix Sq > 0 such that 1 — eo||^/iO||L°° ^ (0) I); ^ ^ [0)^o]- On the 
one hand, for f = g, we deduce that n = v and that uniquely dehnes the 
mapping (pe[g] '■= On the other hand, the function is Lipschitz continuous 
because of □ 

We also recall the following classical Gronwall’s type lemma. 

Lemma 2.9. Assume that u G ^([O, cxd);M+) satisfies the integral inequality 

u(t) <Cie°'^ uq +C 2 f f ds, Vt > 0, 

Jo 

for some constants Ci > 1, C 2 ,uo > 0 and a < 0. Under the smallness 
assumption 

a + 2 G 2 no < 0, 


there holds 

Proof of Lemma VIM We fix A G (OiMO) 2GiMo); so that Ciuit) < A at least 
on a small interval, that is for any f G [0,r], r > 0 small enough, and then 
the integral inequality implies on the same interval 


u{t) < Mo + C 2 Cf^A / ds 


The classical Gronwall’s lemma (for linear integral inequality) and the small¬ 
ness assumption a -|- C 2 Cf^A < 0 imply 


u{t) < Cl Mo 'Ti < Gi Mo < A 


on that interval. By a continuity argument, the first above inequality holds 
on M+ and then with A := CiUq. Next, replacing that hrst estimate in the 
integral inequality we started with, we get 

u(t) < Mo + C2C^uy^ [ Wt> 0, 

Jo 

from which we immediately conclude. □ 

We come to the proof of our main result Theorem ll.2l in the case without 
delay. 
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Proof of Theorem li.i^ in the case without delay. We split the proof into two 
steps. 

Step 1. New formulation. We start giving a new formulation of the solutions 
to the evolution and stationary equations in the weak connectivity regime 
e G (Oj^o], where Sq is dehned in Lemma I2.81 For a given initial datum 
0 < /o G L^(M+) with unit mass the solution / G ^([O, oo); L^(M_|_)) to the 
evolution equation ffTTr]) and the solution to the stationary equation ( [TTgj ) 
clearly satisfy 

dtf + d^f + a{e(p[f])f = 0, /(f, 0) = ip[f{t, •)], 

d^F + a{eM)F = 0, F{0) = M = ^[F], 

where here and below the £ and x dependency is often removed without risk 
of misleading. 

We introduce the variation function g := f — F which satishes the PDF 
dtg = -d^g - a{eM)g - ea\eM)FM[g] - Q[g] ( 2 . 36 ) 

with 


Q[g] := a{e(p[f])f - a{e(p[F])F - a{e(p[F])g - ea'{e(p[F])FM[g], 

where Ai = Aie is dehned in ( p7^ . The above PDE is complemented with 
the boundary condition 

= (p[f{h-)]-(p[F], 


and we may write again 

poo poo 


= j [a{eM)g + ea'{eM)F Ai[g]) + j Q[g]dx 
= Ai[g\ + Q[g], Q[g\ := {Q[g]). 


As a consequence, we have proved that the variation function g satishes the 
equation 

dtg = Aeg + Z[g], Z[g] := -Q[g\ + 5tiQ[g]. ( 2 . 37 ) 

Step 2. The nonlinear term. On the one hand, we obviously have 

{Z[g])=t), \fgeM\^+). ( 2 . 38 ) 
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On the other hand, in order to get an estimate on the nonlinear term Z[g], 
we introduce the notation 


tlj{u) = a{x,emu)fu, 

where, for some fixed g G P(M+), (g) = 0, we have set 

f-=F + g, fu := uf + {1 - u)F, := 

We first notice that ■0(0) = a{eip[F])F and ■0(1) = a{eip[f])f. Second, we 
have 

+ ae{mu)g. ( 2 . 39 ) 

In order to compute we differentiate with respect to u the identity 


mu= aeirriujfudx, 


and we have 


m,, = 


a'^{mu)fndxm'^+ / ae{mu) gdx, 


which implies 


= (1 - / a'^{mu)fudx) ^ / ae{mu) gdx. ( 2 . 40 ) 


We may thus observe that uIq = M.[g], so that '^'(O) = a'^{M)FMe[g] + 
as{M)g, and therefore 

Qia] = ' 0 ( 1 ) -' 0 ( 0 ) -' 000 ). 


Third, from ( [ 2 . 39 ! ), we have 

^/>"( m ) = a"{mu)fu + 2a'^{niu)gm'^ + a'^{mu)fum'^, 

and from (| 2 . 4 oj), we have 
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In the small connectivity regime £ G (0,£o]) ^o||o1loo < 1, we get the bound 


\W{u)\\x < ||a"||oo + 2||a'||oo||fi'||x|m'„| + II 4 II 00 

|2 ll^ll 


< e" 


(1 


+2e^ 


- £||a'||oo)^ 

lla'iP Hall 

11^ llooll^lloo 

(1 - £||a'||oo)^ 


2 


X 


+ 2^1 


00 00 


X 


+ 2 £'^ 


£||a'||oc 
11 ^ ||oo||^ ||oo||^||c 

(1 - ella'lloo)^ 


< eK\ 

for some constant K G (0, cxd). Using the Taylor expansion 

<5[£/] = ^(1) - V’(O) - V’'(0) = f {I-u)ij"{u)du, 

Jo 

we then obtain 




2 

X 


\\Z[g]\\x < 2\\Q[g]\\x < [\l - u)\W\u)\\xdu < C ||^?||^. 

Jo 

Step 3. Decay estimate. Thanks to the Duhamel formula, the solution g to 
the evolution equation (g^Tj) satishes 

9(t) = S^Mifo -F)+ f S^^it - s)Z[g{s)] ds. 

Jo 

Using Theorem 12.71 and the second step, we deduce 

llsWIU < Ce‘>‘||9olLv+ rCe"(‘-*l||Z|9(i>)]|Uds 

Jo 

< Ce-^\\go\\x + CeK f \\g{s)\\lds, 

Jo 

for any f > 0 and for some constant (7 > 1, a < 0, independent of £ G (0, £o]- 
Observing that ||fi'(t)||x G ^([O, 00 ), we conclude thanks to Lemma□ 


3 Case with delay 

This section is devoted to the proof of our main result. Theorem 11.21 in the 
case with delay by following the same strategy as in the case without delay 
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but adapting the functional framework. The main difference comes from the 
boundary term and it will be explained in the first subsection. We have 
already proved in Theorem 11.11 the existence of a unique stationary solution 
(Fg, Me) in the weak connectivity regime and we may then focus on the 
evolution equation. 

3.1 Linearized equation and structure of the spectrum 

In order to write as a time autonomous equation the linearized equation 
fii.i 2 j) - (|i.i 3 j )-( |i.i 4 D , we introduce the following intermediate evolution equa¬ 
tion on a function v = v{t, y) 

dtV + dyV = 0, v{t,0) = q{t), v{0,y) = 0, ( 3 . 1 ) 

where y > 0 represents the local time for the network activity. That last 
equation can be solved with the characteristics method 


v{t,y) = q{t - y)lo<y<t. 


Therefore, equation ( ji.iqj ) on the variation n{t) of network activity writes 

poo 

n{t) = V[v{t)], V[v\:= j v{y)h{dy), 

Jo 

and then equation (jiT^ on the variation q{t) of discharging neurons writes 

Qit) = Oe[9{t),v{t)], 


with 

0^[g,v] :=J\fe[g] + K,eT^[v], 

poo poo 

J\fe[g]-= / aeiMe)gdx, Ke := / a'^{Me)F^ dx. 

Jo Jo 

All together, we may rewrite the linear system fii.i 2 i) -( iTTT^ -( iTTT^ , as the 
autonomous system 

dt{g,v) = J^’,{g,v), ( 3 . 2 ) 

where the operator is defined by 





-d^g - a^g 


<F,V[v] 


-dyV 
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with domain 


D{^,) := {(^ 7 ,n) G x a;); . 7 ( 0 ) = n( 0 ) = 0,[gM}. 

where uj{x) := with 5 > 0 defining in ( [TT^ . The associated semigroup 
Sit) acts on 


X := Xi X X 2 := L\R+) x L^(M+,a;). 


Considering the boundary condition as a source term, we also introduce the 
semigroup Sa^ (t) acting on 


X X 2 := M\R+) X M\R+,u) 


X ; 


with the generator = {A\, Af) given by 

^1(9, v) := -0x9 - as9 - + 5x=oOe[g, u], 

^l{9,v) := -dyV+ 5y=QOe[g,v]. 


In a similar way as in section |221 we have 5'yi^|x = 

As a first step, we establish that the semigroup Sa^ has a nice decompo¬ 
sition structure with finite dimensional principal modes and a fast decaying 
remainder term. 

Theorem 3.1. Assume fli.6l) -( iTT7j )- fli.8p and ( [TTgj ). The conclusions of The- 
orem \2.!^ holds true with a** := max{a*, —5} < 0. 

The result is obtained as a consequence of the Spectral Mapping and 
Weyl’s theory developed in [TTl [TOl [TT] and taken over in section 12.21 For 
that purpose, we ontroduce the convenient splitting of the operator A^ on A 
as A^ = + Be defined by 



and 



It is only necessary to establish the following adequate properties of the 
operators and We skip the rest of the proof and refer to the proof of 
Theorem 12.21 for more details. 
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Lemma 3.2. Assume that a satisfies conditions fli.6p - fli.8p and that b sat¬ 
isfies (prgj). For any e: > 0, the operators Ae and satisfy : 


(%) Ae e X 

(ii) Sis^{t) is of -hypodissipative in both X and X; 
(Hi) the family of operators * AeSjs^ satisfies 


IKS'^g * AeSBe){t)\\^{x,Y) < C'e"*, Va > a\ 


for some constant Ca > 0 and with Y := Y 1 XY 2 , where Yi = fl 

L\{R+) andY2 = BV{W+,uj) n L\{R+,uj). 

Proof, (i) It is an immediate consequence of the fact that V G cn), M) 

and J\fe G M) because of fli.Sp and ( [ngj ). 

(ii) Since S'^i is nothing but the semigroup dehned in fi 2 .i 2 p which 
is a*-dissipative thanks to Lemma [23]- (h), we just have to prove the dissipa- 
tivity of the translation semigroup SJ 32 which is dehned through the explicit 
formula [Si 32 {t)v]{y) = v{y — t)ly_t>o- That follows from 


for any n G X 2 and any f > 0. 
(hi) For {g,v) G X, we have 



-T^S'B(f)( 5 (,n)(a;) = '^{x)D{t) + 5:,=oN{f), 
■A^SB{t){g,v){y) = 5y=oKD{t) + 5y=oN{t), 


( 3 - 3 ) 

( 34 ) 


with 


7 ( 0 ;) := k 5^= o - a'{x)Ffix), 



where here and below the £ dependency is removed without risk of ambiguity. 
We observe that 


|/V(«)| < Cai^f“\\gU, 


|D(«)| < Ce“|K’ll*, 
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for any t > 0. We then compute 

poo 

N'{t) = / - t)l^-t>odx 

Jo 

poo 

= / [a’{x) - - t)l^_t>o dx, 

Jo 

poo 

D'{t) = / v{y-t)ly-t>oh'{(ly) = h'*v{t), 

Jo 

from what we get the estimates 

\N'{t)\ < 

\D’{t)\ < e-^*||n|U,. 


Denoting 


we compute 


Ti{t){g,v){x) := {Sbi * Si3){t){g,v){x), 

T2{t){g,v){y) := {8^2 * A^Si3){t){g,v){y), 


Ti{t){g,v){x) = / Sj3i{s){'y{x)D{t - s) + S^=oN{t - s))ds 

Jo 

= f e^'^'^~"'^~A^\'y{x-s)D{t-s) + 5^_s=oN{t-s)l^_s>o)(is 
Jo 

= * b^){x) + e-A-)N{t - 

T2{t){g,v){y) = [ SB2{s)Sy=o{KD{t - s) + N{t - s))ds 

Jo 

= / 6 y-s=o{KD{t - s) + N{t - s))ds 

Jo 

= [KD{t -y) + N{t - y)) ly<t, 

where we use the notation z/ := 7 We next differentiate the above identity, 
and we get 

dxTi{t){g,v){x) = -a{x)e~'^^^\{v * bt){x) + N{t - x)l^<t) 

-e-^i-){j,^b[{x) + N'{t-x)l^<t) 

_e-^(-)zy(a; - t)D( 0 )l,_t>o - e-^(")iV(0)4=i 
+e-^("V(a;)D(t), 

dyT2{t){g,v){y) = -{kD'+ N'){t-y)ly<t-{KD + Nm5y=t 
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All together, we deduce 


||9.r,(()(9.«)Wlk, < e““||(j,t>)IU. 
\\d,Tm9,v){y)\\x, < e“'‘||(9.0ll;v. 

and the similar estimate for \\{Sj 3 * ASj 3 )(t){g,v)\\x- As a consequence, the 
announced estimate holds for the family of operators Sjs * ASb- □ 


3.2 The vanishing connectivity regime 

When the network connectivity parameter vanishes, £ = 0, the linearized 
operator simplihes as 


A [A ^ (-9x9 - a{.x, 0)g + 6x=oOo[g, v]\ 
°\vj -dyV + 6y=oOo[g,v] ) 


(3-5) 


where C>o[ 5 'w] = a(x, 0)g(x)dx. The associated semigroup is 

exponentially stable as shown in the following theorem. 

Theorem 3.3. There exist some constants a < 0 and C > 0 such that 
S(Ao) n Aq = { 0 } and for any {go, vq) G X, {go) = 0 , there holds 


\\SAo{t){9o,Vo)\\x<Ce°‘^\\{go,Vo)\\x, Vf>0. ( 3 . 6 ) 


Proof of Theorem \d.,‘A Since Aq = Aq, from Theorem 12.41 we have already 
proved that g{t) := Sxi{t)go satishes A lAe®* ||(7 o||a’i for any f > 0 

and any a G (a*,0). We then focus on Aq. The Duhamel formula associated 
to the equation dtv = Al{g,v) writes 

^(t) = SB^^{t)vo + [ Ss^it - s)Al{g{s),v{s)) ds. 

Jo 

Using the already known estimate on g{t), we deduce 

= \\v{t)\\x2 < \\SB^^{t)vo\\x2 + [ - s)<^oA/'o[c/(s)]|U2ds 

Jo 

< e-^^\\vo\\x2 +I e-^^^-9Cen\go\Uds 

< Ce"*||(5fo, 'i^o)IU 

for max{a, —5} < a < 0, which yields our conclusion. □ 
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3.3 Weak connectivity regime - exponential stability of 
the linearized equation 

In this part, we shall discuss the geometry structure of the spectrum of the 
linearized time elapsed equation in weak connectivity regime taking delay 
into account and using again a perturbation argument. 

Theorem 3.4. There exists some constants £o > 0, C > 1 and a < 0 such 
that for any e G [0, £o] there holds ^(klg) n Aq, = {0} and 

< C'e"*||(^o,'yo)IU, (3-7) 


for any {go,vo) G X such that {go) = 0. 

We start presenting a technical result needed in the proof below. 

Lemma 3.5. The operator is continuous with respect to e, and more 
precisely 

W^e - Ao\\ss(x) < 0{e). ( 3 . 8 ) 

Proof. For all {g,v) G X, we have 


Ae 



Ao 



f-dxg - a^g - a{FeVs[v] + 6^=oOe[g, 
-dyv + 6y=oOe[g,v] 



f-dxg - a{x, 0)5( + 5 x=oOo[g, 

\ -dyV X 5 y=oOo[g,v] 



( 3 - 9 a) 

(3-9b) 


From ( i 3 . 9 ai )-( | 3 . 9 bD , we deduce 




(a(a;, 0) 


ae)g 


- a'^F^V^lv] + Sx=oiOe[g, v] - Oo[g, n]) 
dy=o{Oe[g,v] - Oo[g,v\) 


We then compute 


\\{A,- Ao){g,v)\\x = ||(a(x,0) - a,)g\\xi + \\a{F,V,[v]\\x, + 2\Oe[g,v] -Oo[g,v]\ 

< 3||(a£ - ao)fi'|Ui + 2\\a{FsVe[v]\\x2 

< 3e:||a'||oo||5'IUi + 2£ai||a'||oo(l - e:||a'||oo)||T£||;til|w|U2 
= Ce\\{g,v)\\x, 


which is nothing but ([ 3 . 8 ^. 


□ 
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Proof of Theorem \3.4\ With the help of Lemma 13.51 we may proceed simi¬ 
larly as in the proof of Theorem 12.71 (see also again P?[ [T5] ) and we conclude 
that 

S(T,)nA« = fe}, 

with < 0{e) and is algebraically simple. We observe that 

= ( d^ip - ae(p + ae{(p{0) + ijj{0)) \ 

^ [fj J ydyf! + Keb'ijj{0) + Kebip{0) - bj a'^Fe (fi dx J ’ 

from which we deduce that yl*(l,0) = 0. Then 0 G S(yl*) and = 0. 
Moreover, the orthogonality condition {go) = {{go,Vo), (1, 0));i:’,a’' = 0 implies 
that the exponential estimate (j^jTTj) holds. □ 


3.4 Weak connectivity regime - nonlinear exponential 
stability 

We finally come back on the nonlinear problem and we present the proof of 
the second part of our main result for the case with delay. 

Proof of Theorem \1.^A in case with delay. We write the system as 

dtf = -da^f - a^{V[u])f + 6o^[f, V[u]] 
dtu = -dyUP 5o0^[f,V[u]], 

with 

,^[/, m] = J a{m)f, T>[m] = J bu. 

We recall that the steady state {F,U), U := Mly>o, satisfies 

0 = -d^F - a,{M)F + 6oM 
0 = -dyU + SoM, M = V[U] = ^[F,V[U]]. 

We introduce the variation g := f — F and v = u — U. The equation on g is 

dtg = -d,g-a,{V[u])f + a,{M)F + 6omfMy^]]-nF,V[U]]) 

= -dxQ - ae{M)f - a'^FV[v] - Q[g, v] + 5oO[g, v] -f 5oQ[g, v] 
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where 


Q[9,v] ■= {Q[9,v]), 

2^[9,v] := -Q[9,v] + SoQ[9,v], 

with Q[g,v] denoting that 

Q[g,v] := a,{M)F-a,{V[u])f + a,{M)f + a',FV[v] 

= $(0)-<h(l)+ $'(0), 

where 

$(0) = a,{v[eu + (1 - e)u]){e / + (i - e)F). 

The equation on v is 

dtv = -dyV + 6o{^[f,V['^]]-^[FMU]]) 

= -dyv+ SoO[g,v]+6oQ[g,v] 

= ^l{9,v) + 2:^[9 ,v], 

where 

Z^[g,v] := 6oQ[g,v]. 

We then write the associated Duhaniel formula 

{9{t),v{t)) = SaM{9o,vo) + [ SA^{t - s)Z[g{s),v{s)]ds. 

Jo 

Refering to Step 2 in the proof of Theorem 11.21 gives a similar estimate 
ll^[ 5 'W]IU < C\\{g, v)\W, we then conclude the rest part of Theorem 11.21 
in the case with delay. □ 
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